arXiv:hep-th/9507059vl 10Jull995 


EXPLICIT RELATIVISTIC VORTEX SOLUTIONS 


FOR COOL TWO-CONSTITUENT SUPERFLUID DYNAMICS 

Brandon Carter and David Langlois 
Departement d’Astrophysique Relativiste et de Cosmologie, C.N.R.S., 

Observatoire de Paris, 92195 Meudon, France. 

Racah Institute of Physics, The Hebrew University, 

Givat Ram, 91904, Israel. 

(February 1, 2008) 


Abstract 

We give a class of explicit solutions for the stationary and cylindrically sym¬ 
metric vortex configurations for a “cool” two-component superfluid (i.e. su¬ 
perfluid with an ideal gas of phonons). Each solution is characterized only by 
a set of (true) constants of integration. We then compute the effective asymp¬ 
totic contribution of the vortex to the stress energy tensor by comparison with 
a uniform reference state without vortex. 

I. INTRODUCTION. 

The subject of investigation in the present work is the class of vortex type (stationary 
cylindrically symmetric) configurations for the relativistic generalisation of Landau’s two 
constituent superfluid theory, using the recently derived Lagrangian formulation in which 
the independent variables are the superfluid phase scalar ip and the entropy current vector 
s^. More particularly, it will be shown that the complete set of such vortex solutions is 
obtainable in analytic form in case of the “cool” phonon dominated limit regime. The 
required form of the Lagrangian for this “cool” regime has been found p| to be simply given 
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by 


C = P-3'iP, ( 1 . 1 ) 

where P is the “cold” pressure function, depending just on the superfluid momentum cov¬ 
ector 


fip Ti^p^p , 


( 1 . 2 ) 


which governs the zero temperature limit for which the entropy current vanishes, while the 
thermal contribution -0 represents the generalised pressure of the phonon gas, which is given 
by a simple algebraic expression (4.3) involving as well as fip and the “cold” sound speed 
Cg determined by the pressure function P. 

It is to be remarked that whereas accurate treatment of non-st at ionary conhgurations 
with non-zero temperature would require allowance for viscosity of the entropy current, 
however, no loss of accuracy will be entailed by the use of the strictly conservative treatment 
as in Landau’s original model for treating the stationary equilibrium conhgurations under 
consideration here. Indeed under these circumstances there will in any case be no dissipation, 
the only effect of viscosity being to ensure that the “normal” part of the how is constrained 
to have a conhguration that is rigid. 

Unlike the neglect of dissipation that is implicit in our use of a strictly Lagrangian for¬ 
mulation, the limitation that the states under consideration here be restricted to the “cool” 
regime governed by a Lagrangian of the particular form (|1 . 1|) in the sense explained above 
represents a signihcant physical limitation. It is provisionally necessary to postpone detailed 
quantitative analysis of states in the physically very interesting “warm” regime nearer to 
the phase transition (beyond which lies the regime of “hot” states in which superhuidity 
is absent altogether). The reason is not merely the expectation that the equation of state 
governing the explicit form of the dynamical equations in the “warm” regime would be too 
complicated to be easily tractible. A more compelling obstacle is the consideration that 


an appropriate “warm” generalisation of the “cool” equation of state (|1.1|) is not not yet 
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available at all. Although extrapolation beyond the “cool” regime is not yet possible for all 
the explicit quantitative results obtained below, it will nevertheless be found possible to 
obtain many useful qualitative results that remain valid even in the “warm” regime, since as 
shown in Section B the relevant conservation laws provide sufficiently many first integrals 
to solve the differential part of the vortex problem completely in the general case. One is 
thus left with a purely algebraic problem that remains intractible in the general case, but 
that is easily solved in the “cool” limit. 

The formalism that will be used in the present work is fully covariant in the general, not 
just special, relativistic sense, but since the kinds of (neutron star or laboratory Helium) 
vortex that are envisaged are of limited scale it will be justifiable to ignore gravitation, i.e. 
to use a simple Minkowski space background, in the actual application to explicit solutions. 
As well as being obviously important for applications such as neutron star interiors 
in which relativistic effects are actually quite large, the use of a covariant formulation of 
superfluidity theory is also advantageous as a source of supplementary physical insight even 
for the familiar laboratory example of liquid Helium-4 in which relativistic corrections are 
quantitatively negligible |^. 

II. THE EQUATIONS OF MOTION FOR TWO-CONSTITUENT SUPERFLUID 

DYNAMICS. 

Quite generally (not just in the “cool” regime) the dynamical equations for the generalised 
Landau theory are determined by the Lagrangian C in terms of the entropy vector and 
the momentum covector Pp given by (|1.1| ) together with their dynamical conjugates, namely 
the thermal 4-momentum covector 0p and the particle number current vector that are 
constructed |l|] from the Lagrangian according to the infinitesimal variation formula 

dC = Qpds^ — n^dup . (2.1) 

In terms of these quantities, the equations of motion consist just of the thermal momentum 
evolution equation 


3 


s'’V[p0<j] = 0 


(2.2) 


(using square brackets to denote index antisymmetrisation) together with the usual particle 
and entropy conservation laws 

Wy = 0 (2.3) 

and 


vy = 0 . 


(2.4) 


The formulation that has just been summarised has the technical advantage of being 
particularly economical in so much as it involves only 5 independent component variables, 
namely the phase scalar ip and the 4 independent components of the entropy current vector 
. This feature of economy has recently been exploited for the purpose of setting up a 
correspondingly economical Hamiltonian formulation of the theory [Q, and it has also been 
exploited as a guide to the formulation of an analagously economical theory for describing 
thermal effects in superconducting cosmic strings 0. 

The dynamic equations (|1.2|) , (|2.2| ), (^^) , (|2.4|) entail the corresponding pseudo (in flat 
space strict) energy momentum conservation law 


VpT^‘" = 0 


(2.5) 


where the stress momentum energy density tensor is given in terms of the corresponding 
generalised pressure function 


T = £-Qy 


( 2 . 6 ) 


by 


n = • (2.7) 

An obviously convenient way of choosing the 3 independent variables in a fundamental 
state function of the form characterised by ( [^.1| ) will be to take them to consist of the thermal 
rest frame entropy density s as given by 
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(2.8) 


= —SnS^ 


the new cross product variable y given by 


cV = 


(2,9) 


together with the effective mass variable y given by 


c^/i^ = 


( 2 . 10 ) 


It can then be seen that the secondary variables and 0p will be given in terms of the 
primary variables Hp and of this formulation by 


n‘ 


' = - Asi") , + ylfip) , 


( 2 . 11 ) 


where the relevant dilation, determinant and anomally coefficients <h^, /C, and A are given 
by the partial differentiation formulae 




IdC 




IdC 


= 


/i 5/i ’ sds ' 

In terms of these variables, (2.7) can be rewritten as 

+ JCs^s^) + . 


2y dy 


( 2 . 12 ) 


(2.13) 


III. THE CLASS OF STRONG EQUILIBRIUM STATES OF A CYLINDRICAL 

VORTEX. 

The class of vortex conhgurations to be dealt with here is of the maximally symmetric 
type, as characterised by both stationarity and cylindrical symmetry. This means that there 
are three independent commuting Killing vector symmetry generators, k^, and say, 
which may be taken to be the generators of time translations, longitudinal space translations 
(parallel to the axis) and axial rotations, respectively corresponding to ignoreable coordinates 
t, z, and (f) say, where the latter (like the superfluid phase variable ip) is periodic with period 
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2ti. Within this stationary cylindrical category, the class of solutions that will be dealt 
with here is characterised by the condition of equilibrium in the strong sense, which is to 
be understood here as implying exclusion of the (experimentally possible, but under natural 
conditions unlikely) presence of any source or sink of energy or fluid flux in the central core 
of the vortex (where the superfluidity conditions inevitably breakdown). 

In view of the conservation laws (2.3) and (2.4), the condition of equilibrium in this 
strong sense evidently entails that the flow has no radial components and so is necessarily 
helical : this means that the current vectors will be confined to the timelike hypersurface 
generated by the Killing vectors, so that they will be expressible in the form 

= i>{k^ + + ojm^) , (3.1) 

= a(F + V^P + VLmP) . (3.2) 

The class of such helical flow conhgurations includes the specially simple and important 
subclass of circular flow configurations, namely those for which the longitudinal translation 
velocities v and V can be taken to be zero (for a suitably adjusted choice of and i^), as 
is the case in the kinds of vortex that occurs most commonly in practice. Whereas both the 
amplitude factors, u and a, and in the particle current case also the translation velocity v (if 
relevant) and the angular velocity u, can be expected to be radially dependent (though they 
must of course be independent of the ignoreable coordinates, t, z, 0), it is however implicit in 
the “strong” equilibrium condition that in so far as the entropy constituent is concerned the 
corresponding “normal” translation velocity V (if relevant), and the corresponding “normal” 
angular velocity hi have to be strictly uniform, i.e. constant in the radial direction as well. 
This latter requirement expresses the condition that the “normal” constituent flux should 
be rigid in the sense of having its current vector aligned with a Killing vector held, k^ say, 
whose explicit form will be given by 

kP = kP + ViP + QmP = a-^sP . (3.3) 

The reason why it is implicit in the “strong” equilibrium condition that the “normal” con¬ 
stituent of the how should be rigid, is that in a realistic treatment “normality” implies the 
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presence of at least a small amount of viscosity which for non-rigid motion would cause 
dissipation. This could be compatible with stationarity only in the presence of a radial flux 
of entropy and hence also of energy, which would need to be supplied by the core even if the 
only heat sink were at the outer boundary. 

As well as the simplihcations involved in (|3.1| ) and more particularly in ( ^.3| ), it is to 
be remarked that the restriction that we are only considering states of equilibrium in the 
strong sense - requiring a non dissipative flow conhguration - has the practical advantage 
that there is no loss of realism in our use of a strictly conservative theory, of the kind 
summarised in Section 1: there would be no point in including the “normal” viscosity terms 
that are relevant in other contexts, because in the states treated here their presence would 
have no effect. 

The condition that the “normal” part of the flow is rigid makes it subject to a relativistic 
generalisation [Q of a (Jacobi type) variant of the classical Bernoulli theorem. The condition 
that a Killing vector, say, generates a symmetry of flow means that for any physical well 
dehned held (though not of course for a gauge dependent quantity such as the superhuid 
phase variable ip) the corresponding Lie derivative should vanish. In the case of the thermal 
momentum covector 0p this symmetry condition is expressible as 


2FV[p0.] + V.(F0p) = 0 . 


(3.4) 


The condition that the how satishes a rigidity condition of the form ( p.3|) can now be used 
to rewrite this Lie invariance condition in terms of the relevant current vector in the more 
specialised form 


2s^V[p0.] = -aV.(F0p) . (3.5) 

Up to this point the reasonning has been purely kinematic. However if we now use the fact 
that the momentum covector in question is subject to a dynamical equation of the standard 
form ( |2.2|) it can be seen that the left hand side of ( ^.51) will simply vanish, leaving us with 
a uniformity condition that provides us with a hrst integral of the motion in the form 
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F0 =-0, 


(3.6) 


where 0 is a Jacobi-Bernouilli type constant, which is interpretable as an effective tem¬ 
perature that (like the qnantities V and introdnced above) is uniform thronghont. The 
existence of this nniform effective temperatnre (as measnres with respect to the rigidly coro¬ 
tating frame of the “normal” flow) is interpretable as the relevant application of the “zeroth 
law of thermodynamics”. 

Althongh the particle cnrrent does not satisfy a rigidity condition, the fact that the 
associated dynamical eqnation is not jnst of the standard form bnt has the more restrictive 
form of the irrotationality condition ( |1.2| ) provides ns with not merely one bnt three more 
independent Bernonilli type constants, corresponding respectively to the three independent 
Killing vectors. One way to see this is to start by snbstitnting pp in the place of 0p in 
the Lie invariance condition (|3.4|) . Since it is obvious that the dynamic condition (1.2) 
will then immediately anullate the hrst term, one is again left with a uniformity condition 
that provides another Bernonilli type constant, E, which is interpretable as the energy per 
particle with respect to the rigidly rotating “normal” frame, and which is given by 


fc^/ip = —E. 


(3.7) 


Furthermore, since this latter result is not of the Jacobi type in that its derivation did 
not actually depend on the rigid corotation property that singles out the combination 
it is clear that an analogous, ordinary rather than Jacobi type, Bernonilli constant can be 
obtained by replacing with any other symmetry generating Killing vector. We thus obtain 
not merely one but three more independent Bernonilli type constants, E, L, and M say, 
corresponding to the three independent Killing vectors, £^, and m^, and interpretable 
respectively as representing the energy, longitudinal momentum, and angular momentum, 
per particle, which will be given by 


fc^/ip = —E , = L , m^j^p = M . 


(3.8) 


The prototype example (|3.7|) is of course not independent of these, but is given in terms of 
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them as the constant coefficient linear combination 


E = E -VL-QM 


(3,9) 


Up to this point everything that has been done is fnlly general relativistic in the sense of 
being applicable with respect to an arbitrary system of coordinates {x^} for an arbitrarily 
curved space time metric . If we wished to allow for self-gravitation we would still be 
left with a non-trivial differential system of Einstein type equations to be solved. However 
if, as will be sufficient in all the most obvious applications, we are willing to treat the 
gravitational spacetime background as hxed in advance, independently of any feedback from 
the superfluid system, then it is apparent that the work carried out already above has been 
sufficient to fully integrate all the relevant differential equations. To see this, it suffices to 
notice that, as the helicality property implies the absence of any radial components, each 
of the pair of currents involved has only three instead of four independent components, 
which means that the system is characterised locally by a total of only 6 independent held 
components. Therefore, in order to solve the corresponding differential system, it suffices to 
obtain a corresponding set of just 6 hrst integral constants. Precisely such a set is provided 
by the work that has just been described: assembling them all together the complete sextet 
of constants can be listed as the temperature 0 and the linear and angular velocities V 
and H associated with the “normal” part of the system, together with the energy E (or its 
corotating counterpart E) and the linear and angular momenta L and M associated with 
the superhuid part. 

It is to be remarked that as far as the physical characterisation of the solution is concerned 
only 5 of these constants are independent, since by a longitudinal Lorentz adjustment either 
V or L given any desired value without loss of generality, the most mathematically convenient 
choice in the present formulation being to use the “superhuid frame” in which L = 0. 

It is also to be remarked that whereas the hrst hve members of the sextet that has just 
been listed can take values in a continuous range, the quantisation condition on the phase 


variable (f in (|1.2|) entails of course that the last one, M, should be restricted to integral 
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multiples of ft. Moreover, since vortices with higher values of the relevant winding number 
will typically be unstable, it is in common circumstances sufficient to consider only the 
almost trivial case with M = 0 (which might be deemed unworthy of description as proper 
vortex) and the hrst proper vortex possibility as given (subject to a judicious choice of sign 
convention) hj M = h. 

Although, in a given background, the differential part of the problem is in principle 
solved by the foregoing derivation of the sextet of constants 0, V, ff, E, L, M, there 
remains what - depending on the complexity of the relevant equation of state - may be a 
highly non-trivial algebraic problem that still has to be solved if one wishes to obtain all the 
other important non-uniformly distributed physical held quantities involved in explicit form 
as variable functions of a suitably dehned cylindrical radial coordinate r say, which can be 
taken to be the only independent variable in each solution. 

Since such superhuid vortex conhgurations as can in practice be set up artihcially under 
laboratory conditions, or as can be expected to occnr naturally in neutron stars, are all 
characterised by dimensions that are small compared with the length scales over which self 
gravitational effects become important, they should be adequately describable in terms of a 
local gravitational background that is not merely given in advance (which is all that is needed 
for complete integrability as described above) but that can be taken more particularly to be 
flat. This allows us to work with a system of the standard cylindrical form {x°, = 

{t, z, (j), r}, snch that only the radius coordinate r is non-ignoreable, while the others are the 
ignoreable coordinates t, z, 0, mentionned above in association with the Killing vectors k^, 
E, so that the metric will be given by the familiar formula 

gpcrdx^dx"^ = —c^dt^ + dz^ + r'^df^ + . (3.10) 

With respect to this coordinate system the independent Killing vector symmetry generators 
will be given by 

^ (1,0,0,0) , (0,1,0,0), ^ (0,0,1,0) . (3.11) 
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The corresponding components of the primary momentum covector and current vector in 
terms of which the formulation based on the Lagrangian C is formulated will be given by 


/ip {—E, L, M, 0) , cr(l, V, n, 0) . 


(3.12) 


It is thus apparent that, in order to obtain the complete solution for all the relevant dy¬ 
namical variables, all that is still needed is to find the functional dependence of the single 
variable a on r. 

It can immediately be seen that the required function a is directly proportional to the 
quantity as dehned by which can be evaluated using (|3.12|) simply as 



(3.13) 


It can similarly be seen that a is also related by a simple though radially variable propor¬ 


tionality factor to the quantity s as dehned by (|2.8|) , which can be evaluated using (|3.12|) 


as 



(3.14) 


It follows that although it is still necessary to hnd out how a depends on r in order to obtain 


the corresponding radial dependence of the separate scalar state variables s and y, the ratio 


y^!s has a radial dependence that is immediately expressible in terms of the Lorentz factor 


t/2 ri2„2 _i/2 




(3.16) 


associated with the rigid velocity distribution of the “normal” how by the formula 



(3.16) 


s 


Although ( |3.12| ) provides only the ratio (|3.16|) , but not the absolute values, of s and y, it 
is more helpful for the third of the scalar helds that is needed for evaluating the Lagrangian 
function, namely the ehective mass function y as dehned by (|2.1U| ), which can be seen to be 
given explicitly as a function of r by the easily memorable formula 
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(3.17) 


2 2 t2 

It is evident that in the “cold” limit for which the entropy cnrrent vanishes the explicit 
solution is directly available in terms of just the 3 constants involved in (p.l7|) namely E, M 
and L (of which the last is physically redundant, being adjustable to zero by a longitudinal 
Lorentz transformation) since in that case the Lagrangian depends only on the single scalar 
/i. It is noteworthy that while the very simple formula (|3.17|) provides the complete solution 
to the (strong) equilibrium problem for a relativistic vortex in the cold limit, it remains 
formally valid without change, although it is no longer sufficient to provide the complete 
solution all by itself, in the generic “warm superfluid” case. 

To complete the solution of the vortex problem in the “warm” case, we must make use of 
the remaining (sixth) constant of integration, which (unlike the other hve) does not appear 
in (|3.12| ), namely the corotating “normal frame” temperature 0 as introduced in ( |3.6| ). With 
the aid of ( p.ll| ) it can be seen that (|3.6|) can be rewritten in the form 


0cV = Ei^^-C) , 


(3.18) 


where the bracketted function on the right is given explicitly by 


dC ydC 

T — T = —s-- — 

ds 2 dy 


(3.19) 


In principle, since y is given directly by (|3.17|) , all that remains to be done to obtain the 
other two state variables y and s (and hence everything else) as functions of r is to carry 
out the simultaneous solution of the pair of algebraic equations ( |3.16| ) and ( |3.18|) using the 


explicit radial dependence given by ( |3.15| ) and ( p.l7| ). In practice however, whereas the 


solution is immediate in the cold limit (for which y and s vanish so that (|3.17|) suffices by 
itself), in the generic “warm” case it will not even be possible to begin to tackle the problem 
until the form of the state is explicitly known, and even then it might reasonably be feared 
that the ensuing non-linear form of the functional form of the expression given by (|3.19|) 
would make it hopelessly intractible. What will be shown in the next section is that in the 
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“cool” limit as described by the equation of state obtained in the preceeding work 0, the 
particular form of the function given by side of (|3.19|) will be such that the solution will, 
rather surprisingly, be obtainable in an explicit analytic form. 

IV. THE EXPLICIT SOLUTION FOR THE COOL LIMIT CASE. 

In the “cool” limit for which the entropy current is describable simply as a gas of phonons, 
the corresponding limit form of the equation of state that is obtained as the relativistic 
generalisation of the corresponding classical formula of Landau P] is given by an expression 
of the form 


C = P — 3i/j , 


(4.1) 


in which P is the zero temperature pressure function depending only on the single variable 
fi and determining the corresponding zero temperature sound speed Cg by the formula 

-2 d/i (PP 


= h- 


(4.2) 


dP djj? ’ 

and ip is the generalised pressure function of the phonon gas, which is expressible as a funcion 
of all three of the independent scalar variables s, y, and /i by the formula 




(4.3) 


where h is identihable with sufficient accuracy for most purposes with the usual Dirac Plank 
constant h, its exact value being given by 


47r V 6 / ^ ^ 

It can be seen to follow immediately, just from the fact that this function ip is homoge¬ 
neous, of order 4/3 in the variables s and r/^, that the function given by ( p.l8| ) will work out 
in this “cool” limit case simply as 


-C = Aip . 


(4.5) 
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Using (|3.16|) , it can be seen that the equation obtained by substituting ( |4.5|) in (|3.18|) can 
be solved explicitly to give the formula 


.30x3/2 , 


C2 -c2)E2 




) . 


(4.6) 


whose right hand side is interpretable as a function just of the pair of variables F and p, since 
the latter determines by the relation (|4.2|) that is obtained from the zero temperature 
equation of state. The formula (|4.6|) thereby gives y directly as an explicit function of 
the radius r since both F and p are already known by ( p.l5|) and (|3.17|) as functions of r. 
Having thus evaluated y one can then immediately obtain the remaining state variable s as a 
function of r using ( p.l5|) and (|3.16|) , thereby completing the solution of the vortex problem. 

The simplest kind of zero temperature equation of state that can be envisaged for illus¬ 
trating the application of this formula is the polytropic kind which, for a given “rest” mass 
m per particle, takes the form P cx (p — where the index a is a positive constant 

(taking the value a = 3 in the standard case of a relativistic gas with kinetic energy large 
compared with the rest mass energy), so that, by (E3)> one obtains = (c^/a)(l — m/y), 
which conveniently reduces to a constant in the high density limit for which y becomes large 
compared with m. 


V. THE NATURAL CUT OFF RADIUS FOR A VORTEX CELL. 

The validity of the solutions obtained in the previous sections is of course limited to a 
hnite range of the radial coordinate r which cannot exceed the critical null corotation radius 
at which the rigid rotation velocity reaches the speed of light so that the Lorentz factor F 
dehned by (|3.15|) becomes inhnite. There is no danger of approaching such a singularity in 
the usual laboratory experiments on exactly cylindrical superfluid vortices for which wall of 
the container limits the radius to a value that is very small compared with the value. The 
occurrence of such a singularity is also avoided in more natural contexts (such as that of 
the neutron star material whose analysis is the ultimate purpose of the present work) under 
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conditions such that the local superfluid flow can be represented by a honeycomb lattice 
of hexagonal vortex cells. In such an application a cylindrically symmetric solution of the 
kind obtained in the previous sections can be used as approximate description of the motion 
within an individual hexagonal vortex cell. In a stationary state the honeycomb lattice 
will be in a state of rigid rotation with an angular velocity that must be the same as 
that of the normal constituent (if any) in order to avoid the presence of dissipation (which 
would be incompatible with strict stationarity). In such a configuration, each approximately 
cylindrical vortex cell will have a natural outer cut off radius, 

r = f (5.1) 

say, where f is just half the distance between neighbouring vortex cores. (The axial symmetry 
approximation should be extremely accurate in the inner region r « f where most of the 
vortex energy is concentrated, so the errors due to neglect of the breakdown from axial to 
hexagonal symmetry in the outer regions should be relatively small.) The symmetry between 
neigbouring vortices of the lattice implies that at the midpoint between nearest neighbour 
vortex cores the angular velocity, not just of the normal flux s^, but also of the particle 
flux should should agree with that angular velocity of the lattice, with respect to the 
frame in which the longitudinal velocity V of the rigid normal flow vanishes. This means 
that provided we fix the choice of the longitudinal Killing vectors and by imposing the 
requirement 


V = 0 


(5.2) 


then the appropriate cut off radius (|5.1|) can be characterised simply by the condition 


UJ = hi , 


(5.3) 


where u is the radially variable angular velocity of the particle current as introduced in 
(EH). Since it can be seen from ( |2.11|) and ( ^.121) that this angular velocity variable will be 
given by 
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a;r^ M — ilr'^Aa 

A E — AAa 


(5.4) 


the equation (|5.3|) is easily soluble. The natural cut off radius value is thereby found to be 
given by the simple formula 

AM 


A = 


En 


(6.6) 


This formula can be used to relate the angular velocity a; of the rigidly rotating frame to 
the corresponding circumferential value w of of the circulation w per unit area. Since the 
momentum circulation k say will be given in terms of the angular momentum constant M 
hy K = 2 ttM, the corresponding mean vorticity, w for a circuit of radius r will be given by 


w = 


2M 
^2 ’ 


(5.6) 


SO that for the corresponding effective macroscopic vorticity, as given by the momentum 
circulation per unit area for each entire vortex cell, one obtains the formula 


w = 2EVt/A . 


(5.7) 


For the purpose of applications to the analysis of macroscopic properties of a vortex 
lattice, it is of particular interest to evaluate total integrated quantities, such as the total 
effective longitudinal energy momentum tensor of a vortex as functions of the corresponding 


total longitudinal currents, within the relevant cut off radius (|5.1|) . 

It is shown in the appendix how the asymptotic deviation of the sectionally averaged 
value Q of any suitable physical quantity Q from its value in an appropriately chosen 
homogeneous reference state goes like Inr for large cut-off distances r, the multiplicative 
factor depending on a corresponding “net asymptotic deviation coefficient” Q. In order for 
the averaging process to be meaningful, it is not necessary that the physical quantity Q under 
consideration should be a scalar in the 4-dimensional sense, but is sufficient that it should 
be a scalar with respect to the coordinates 0 and r of the section at constant t and z over 
which the integrals are taken, as is the case for tensor components as longitudinally projected 
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onto the symmetry axis whose coordinates are the subset {x*} {t, z}, for {i} = {0,1} 

within the full set {x^} introduced in ( p.lOj ). This applies in particular to the corresponding 
longitudinally projected components of the stress tensor, whose sectionally integrated total 
will be given in terms of the corresponding averaged values T\ by 

(Tt) = ■ (8.8) 


The foregoing formulae are also valid for the trace of the orthogonally projected part of the 
stress tensor, which is twice what is known as the the lateral pressure U say, whose dehnition 
is expressible equivalently by 


277 = - T\ , (5.9) 

so that its average will be expressible in terms of the scalar functions introduced in Section 
1 by 

277 = 3T-A-^. (5.10) 


According to the general formula obtained in the appendix, the averages in 
given asymptotically by an expression of the form 


TJ 


T\ ~ 

ek 


Mw , ftry I 
n-^ln{^|, 

^ 111 J 


will be 


(5.11) 


as w —0 (using the standard notation convention according to which the symbol ~ relates 
quantities whose ratio tends to unity in the limit under consideration) where are the 
corresponding reference state values, and is a hxed vorticity value determined according to 
( b.6|) by a suitably chosen hxed “sheath” radius whose precise specihcation is unimportant 
when w is sufficiently small. 


VI. STRESS - ENERGY COEFFICIENTS FOR A COOL VORTEX CELL. 

The purpose of this section is to use the procedure described in the appendix to evaluate 
the coefficient T^, taking the variables q that hx the reference state to be simply the longitu¬ 
dinal components /i* of the particle momentum (the other components being irrelevant since 
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the axial symmetry ensures that they average to zero ) and the longitudinal components s* 
of the entropy current. This choice for the reference state is not at all compelling, and we 
shall consider for instance in the next section, in the cold limit, a reference state dehned by 
the longitudinal components of the particle current. Our choice here is simply motivated 
by the computational simplicity when one will wish to apply the sectional averaging to a 
vortex solution, for which the most convenient variables are the particle momentum and the 
entropy current as illustrated in Section JTT[ 

The reference state is thus chosen to satisfy the condition that the corresponding reference 
momentum components /Xig and reference current components agree respectively with 
the average momentum components Jti and the current components s* of the vortex, so that 
in the notation of ( |^) one has 

= 0 , = 0 . ( 6 . 1 ) 

According to ( |A23| ) , this means that the required “net” asymptotic deviation coefficients 
T\ will be given in terms of the corresponding “gross” asymptotic deviation coefficients T\ 
by 

A = A - ;< An. (6.2) 

using the fact that 


fjji 0 , 


(6.3) 


which follows obviously from (|3.12|) . 

To apply the formula ( |6.2|) we hrst need to use ( |A19| ) to obtain the required “gross” 
asymptotic deviation coefficients s* for the entropy current itself, whose “net” asymptotic 
deviation coefficients will automatically vanish s* = 0 as an automatic consequence of the 


choice ( |6.1| ). In the case of the vortex solution, the dependence of s* on p is completely 
conhned to the coefficient a as can be seen from ( |3.12|) , which is itself related to the “cross” 
scalar according to ( p.l3| ), so the deviation formula (|A19|) gives 
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gl = 


/i 


y"^ djj, 


(6.4) 


where the limit r —> cx) is implicit in the right hand side, as will be the case in all the 
formulas giving asymptotic coefficients. 

We are now ready to consider the application to the stress tensor which can be written 
in the form 


r', = (4‘V. + + [£ + (sX + y'‘A)] g 


(6.5) 


It is now straightforward to compute T\ by using systematically the dehnition (|A19|) . 
We skip here the details of the calculation, inviting the reader to consult the similar but 
simpler procedure for the cold limit that we will detail in the next section. The result 
involves many terms. It turns out that when one adds to the expression for T\ the extra 
term in order to get the net asymptotic coefficient T\, a lot of cancellations occur, some 
due to the relation 


ds = ^dy^ + ^sTH{T-^) , 
so that the hnal result can be simply expressed as 

, , 2 2 /' .,2 , 2S'S>"'; 2 


( 6 . 6 ) 


(6.7) 


We now consider the asymptotic behaviour of the lateral pressure contribution (|5.10|) 
which can be seen to be given locally by 


n = C + sX + y^A + i(2i^ - - r =) 


( 6 . 8 ) 


Although the expression that one obtains “naively” for 77, i.e. by just applying the dehnition 
( |A16|) , involves a lot of terms, it can be checked, by using (|6.6|) as well as the differential of 
the relation (p.l8|) written in the form 


o2 0 

$2^ + $2^= = , 
7/2 E 


(6.9) 
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whose right hand side is constant, that the “gross” asymptotic deviation coefficient cancels 
out altogether, i.e. one simply gets 

77 = 0. (6.10) 


This result could of course have been predicted in advance from the requirement of overall 
stress balance on the outer boundary. 

Despite of the stress balance condition (|6.10|) , the lateral pressure will nevertheless have 
a non vanishing “net” asymptotic deviation coefficient that will be given by 


77 = /i s 




y 




( 6 . 11 ) 


y ds dy"^ J 

The preceeding results are valid for any equation of state C. We now specialize these 
results to the “cool” equation of state ( [4.1|) for which we have explicit vortex solutions. In 
this case, the dilation coefficient is given by 

1 dP ?, d'lp 


= 


( 6 . 12 ) 


dy yc^ dy ’ 

the second term on the right hand side being homogeneous of order 4/3 in the variables s 
and y"^. It can thus be checked easily that the lateral pressure is simply given by 


77 = 


4 dip 
(? dy 


(6.13) 


whereas the coefficient T\ is given by substituting ( |6.12| ) in (|6.7| ). 


VII. STRESS - ENERGY COEFFICIENTS FOR A COLD VORTEX CELL. 


In this hnal section, we restrict our attention to the case of the cold limit in which there 
is no entropy vector. The superfluid in this zero temperature limit reduces to a particular 
case of perfect fluid and the Lagrangian C reduces simply to the pressure function P{y). In 
this simple case the current will by given, in terms of the “dilatonic” amplitude $ that plays 
a key role in the vorticity variational formulation of the zero temperature fluid model, 
just by 
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nP = , 


(7.1) 


where and its derivative are given by 


(|)^ = 


= 4 (^-l) 


/r ' dfj, 

and where the particle nnmber density n is dehned by 


(7.2) 


nPup = —c^n^ 


(7.3) 


and can also be derived from the pressnre fnnction P(/i) by 

,dP 


n = c 


- 2 ' 


d^ 


(7.4) 


One can also introdnce the energy density of the flnid, p(n), as the Legendre transform of 
the pressnre fnnction P(p), 


p{n) = pn — c ^P{p) 


(7.5) 


It follows from ( p.l2|) that the longitudinal and angular velocities v and uj introduced in 
( p.l|) and the corresponding amplitude factor v will be given by 

c^L c^M E 


V = 


UJ = 


^ = 77 ®^ 


(7.6) 


E ’ " ’ c 

Finally, in this zero temperature limit, the stress tensor has the simple perfect fluid form 


T' = + Pj', . (7.7) 

It is evident from ( |6.2|) that in the cold limit for which the entropy current vanishes 
the “gross” deviations (with respect to the large distance limit) are the same as the “net” 
deviations with respect to a reference state of the kind postulated in the preceeding section. 
This is because the cold limit of such a reference state limit is characterised just by the 
longitudinal momentum components, which are radially uniform, so that deviations from 
their large distance limit values simply vanish; 

= 0 . (7.8) 
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Using (|7.2|) the “gross” variation coefficients for the longitudinal stress energy tensor in the 
cold limit can therefore be obtained simply by setting the entropy current s* to zero in (|6.7|) 
which gives 

T\ = 1^1 - - (pc^ + P)gl (7.9) 

while by (|6.9|) the corresponding “gross” lateral pressure simply vanishes. 

Introducing the longitudinal unit flow vector, say, and the corresponding spacial 
projection tensor 


la = 9^a + 


U^Up = — c^. 


(7.10) 


according to the component prescription g,Q{u°,u^,u'^,u^} = {JP, pT, 0, 0} (in the coordinate 
system introduced in (p.l0|) that we have been using) so that the uniform asymptotic limit 
state is characterised by the condition that its momentum should simply have the form 


the meaning of the formula ([7.9| ) can be interpreted as follows. In terms of the corresponding 
pressure and energy density given, via the equation of state ( [^.1| ) by 

Pe = ^{he} , ^ 

the stress energy tensor of the uniform asymptotic limit state will be given by 


rpp 


PgU'U, + ■ 


(7.13) 


The signihcance of ([7.9|) is that that it specihes the coefficient in the asymptotic formula 

ttM^ 




■% 

w 


(7.14) 


for the difference of the total sectionally integrated stress momentum energy density compo¬ 
nents from the values they would have for the uniform asymptotic limit state (if the vortex 
were absent) as a function of the cut off radius r as given in terms of the corresponding 
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circulation per unit area w by (|5.6|) . The components in this expression are just the 
longitudinal subset of the 4-dimensional set given by ([7.13|) while the formula (|7.9|) for the 
coefficients T\ can be rewritten equivalently, in a form more directly comparable with ( |6.13|) , 


as 


n = - cfyut - ji). 


(7.16) 


Substituting this in (|7.14|) gives the simple explicit formula 

(n) - (ri) ~ ^i>h(2c-^ - cf)uiu, - yj 1.1 {^} 


(7.16) 


The final formula ( |7.16|) is interpretable as meaning that relative to the reference state 
labelled by q, with the same longitudinal momentum component values jUj, the vortex cell 
has an effective energy per unit length, [U] say, that is given by 


PI 


2-3 nn 


(7.17) 


where, [T]^ is the corresponding effective tension, which, for the usual case with the lowest 
quantum value, M = h, will be given as a function of the mean vorticity w defined by (^.61) 
and the amplitude defined by ( |6.2|) as 


[T] - 


.$2 In 1^1 , 
'■ w J 


(7.18) 


wherein - in view of the fact that higher order corrections have been neglected - it is 
sufficiently accurate to take the quantities involved to be those characterising the simple 
vortex free reference state. 

In the non-relativistic limit, in which one can make the substitution 4)^ ~ p/m? where 
m is the relevant fixed Newtonian mass per particle, the expression (|7.18|) can be seen to 
agree with the well known formula given by Hall for the effective vortex tension defined 
as the integrated longitudinal pressure deficit relative to the asymptotic value. Hall pointed 
out that this value is the same as that given by the Feynman formula for the energy per unit 
length, as evaluated in the incompressible case. It is to be remarked that the incompressible 
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Newtonian case corresponds to the non-relativistic limit of the “stiff’ case characterised by 
Cg = c^, and that for a “softer” material the energy per unit length given by (|7.16|) will be 
relatively diminished, so much so that it will actually be negative whenever < c^/2, as 
will typically be the case in realistic applications, including all except perhaps the most 
central regions of neutron stars. This energy dehcit is simply due to the central matter 
dehcit resulting, unless the material is sufficiently stiff, from the centrifugal effect. It will be 
seen in the next section that, when normalised with respect to a hxed amount of conserved 
matter, the energy per unit length will always be positive, and furthermore that its value 
will remain equal to that given by (|7.18|) even for material that is quite “soft”. 


VIII. STRESS - ENERGY WITH RESPECT TO THE CURRENT FLUX 

REFERENCE STATE. 


Rather than the prescription of Section which as we have seen would give no difference 
between “gross” and “net” deviations in the cold limit, let proceed with the analysis of the 
cold limit on the base of a prescription that is less mathematically trivial. Instead of merely 
choosing the reference state to be specihed by the longitudinal components of the momentum 
covector, it is more useful for many physical purposes to choose a reference state to be hxed 
by the longitudinal components of the particle current hux (the other components being 
irrelevant since the axial symmetry ensures that they average to zero). To avoid confusion 
with the previous choice of the reference state that was labelled by the symbol e, we shall 
use a slightly modihed symbol, 0 , to label this new reference state which is chosen so as 
to satisfy the condition that the corresponding reference current components agree with 
average current components n* of the vortex, so that in the notation of (|A8|) one has 


S^n‘ = 0 . 


( 8 , 1 ) 


According to the appendix, in the same way as in the expression (|6.2|) , the relation 
between the “net” and “gross” asymptotic deviation coefficients will be given in this case by 
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( 8 . 2 ) 


T\ = T{ - n* . 

To use this formula we need to the “gross” asymptotic deviation coefficients n* for the 


current itself. Using the last of the expressions (|7.6|) , it can be seen that the asymptotic 
deviation formula ( |A16|) will give 

o2 


rv = 


(1 - . 


(8,3) 


in which the ratio of the light and sound speeds is given by ( |4.2|) . 

In view of the fact that, by ( p.l2| ), the relevant components of the momentum covector 
/ip are constant, as also are those of the metric, they can be taken outside the averaging 
process so as to allow us to write 


T\ = + Pgi 


(8,4) 


Therefore to complete the evaluation of T\ all that remains to be done is to obtain the 
asymptotic deviation coefficient of the pressure scalar P. Applying (|A16| ) and using ([f.4|) 
and (|7.5|) it can be seen that the required result is given simply by 


P = -(pc^ + P) 


(8,5) 


Putting this together with ( |8.3|) in ( p.4|) we obtain a direct derivation of the complete 
expression (|7.9| ) for the “gross” asymptotic deviation coefficient of the stress momentum 
energy density tensor. 

Proceeding towards the evaluation of the “net” deviation coefficients in which we are 
ultimately interested, the next step is to use (|8.3|) again for working out the second term 


required for the application of (|8.2|) . The contribution from the pressure scalar P, which (in 
view of the asymptotic agreement 77 ~ P as /r —/r^) is the same as that for 77, is given 
simply by 




~dP ~dn id 2 


( 8 , 6 ) 


25 















As an immediate corollary we see from ( p.2|) that in to contrast to the case where the 
reference state is specihed by the momentum, the lateral pressure will have a non vanishing 
“net” asymptotic deviation coefficient that will be given by 




(8.7) 


After analogously evaluating the other relevant contributions their combination in the 
complete longitudinally projected stress tensor is obtained as 


Q , ^2 . ^2 


"• a)rn = (^ - + (^ - i)(pc^ + P)<f, 


( 8 , 8 ) 


When this last, rather unwelldy, expression is subtracted off from (O) in accordance with 
the prescription (|8.2|) , we end up with a comparitively simple expression for the required 
“net” deviation coefficient of the stress tensor, which is expressible just by 


n = (i - i)"' 


H ^2 


+ P)g\ , 


(8.9) 


in which, as in preceeding formulae, the relevant values of the quantities involved are those 
of the large distance limit. 

This result is to be interpreted with respect to the reference state as characterised 
by a density given in terms of the average flow by the component prescription 
Mb M®} = {ffo, hJ, 0,0} where is the same unit vector as was introduced in 
(7.10) so that the corresponding current has the form 


K = v” 


( 8 . 10 ) 


In terms of the corresponding density and pressure given, via the equation of state 
(0),by 


P0 = pWI 1 


p _ 


( 8 . 11 ) 


The stress energy tensor of the uniform asymptotic limit state will be given in terms of 
these quantities an expression of the perfect fluid form that is the formal analogue of (|7.13| ) 
namely 
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( 8 . 12 ) 


= p^u'’u, + . 

With respect to this reference state, the expression for the total sectionally integrated stress 
tensor will be given by 

TT A/f‘^ _ / 1 . r m ^ 

( 8 . 13 ) 




This last formula is interpretable as meaning that relative to the reference state labelled 
0 , with the same average longitudinal current vector n®, the vortex cell has an effective 
tension, [T] say, and an effective energy density per unit length, \U] say, that are related 


by 


where \U] is given by 




Pi 


Tih^ 


o fW ^ 

$2ln|^| , 
t w J 


( 8 . 14 ) 


( 8 . 15 ) 


for the lowest quantum value M = h. 

It is to be remarked that relative to the reference state used in the present section, 
as dehned in terms of a hxed quantity of conserved matter, the effective tension can be 
seen from (|8.14D to tend to zero in the limit C 2 « c 2 . It is also to be noted that - as 
a generalisation of the non-relativistic Hall [|^ equality - the last formula (^.151) for the 
energy per unit length is the same as the formula ( [f.l 8 | ) that was obtained for the effective 
tension, albeit relative to a different reference state. However it is only in the “stiff’ case 
characterised by that, as in the incompressible non-relativistic case, the effective 

tension will be equal to the energy per unit length with respect to the same reference state. 


APPENDIX: ASYMPTOTIC FORMS OF SECTIONAL INTEGRALS AND 

AVERAGES. 

The purpose of this appendix is to consider sectional integrals and averages of a generic 
local radially dependent physical quantity Q in the asymptotic limit as the outer cut off 
radius f of the section tends to inhnity, i.e. in the limit as the mean vorticity. 
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w = 


2M 

f2 ’ 


(Al) 


tends to zero. We use angle brackets to denote the sectional integral of Q over the entire two 
dimensional section S say out to the exterior cut off radius f of the cylindrically symmetric 
vortex region under consideration, i.e. we write 

{Q) = J Q dT, = 2 tt J Qrdr . (A2) 

Such an integral will typically diverge in the asymptotic limit f ^ oo, unlike the corre¬ 
sponding average, 

Q = ^ = ^ TQrdr. (A3) 

which will typically tend to a well dehned hnite limit. Provided that Q itself is asymptotically 
well behaved in the sense of tending smoothly to a well dehned asymptotic limit 


(A4) 

as r —*■ cx), its average will obviously converge to the same asymptotic limit i.e. we shall 
have 


Q^Q^- (AS) 

Our interest will therefore focus for the time being on the non trivial difference between the 
average and its limit, as dehned by 

« = (A6) 

and on the corresponding diherence as dehned for the total integral by 

(Q) = (Q) - ■ (A7) 

Whereas it is evident from (|A5|) that 5^Q will tend to zero, on the other hand S^{Q) will 
still diverge, albeit only logarithmicly, i.e. much less strongly than (Q), as f —> oo. 
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In order to analyse this limit it is useful to decompose the range of integration into an 
“asymptotic” zone, f > r > , w < w < and an inner “sheath” zone, r < , w > w^, 

taking the sheath radius to be large compared with the central core radius within which 
the superfluidity breaks down but small compared with the outer cut off radius f. It can be 
seen that the average Q as taken over the entire section out to f will be given in terms of the 
sheath average as taken over just the inner zone extending only to by an expression 
of the form 


Q = —Qe 


w 


"0 ^dw 
V 2 * 


(A8) 


The subject of our ultimate concern is not so much the total integral or average as given 
by the foregoing expressions but rather just the part thereof that is attributable to the 
vortex per se, after subtracting off a contribution that would have been there in its absence. 
What we are specially interested in are not the “gross” differences dehned by (|A4|) but the 
corresponding “net” differences 


S^Q = Q-Q^, (A9) 

where is the value that the quantity under consideration would have in some appro¬ 
priately chosen reference state. We take this reference state not just vortex free but more 
specihcly uniform, which of course entails the exact equality Qq = Q^. In practice we 
shall compute this “net” difference by starting from the “gross” difference and by using the 
relation 


^eQ = (AlO) 

where is the difference between and the aymptotic value in the vortex state, 

i.e. 


=Qe-Q^ 


(All) 


It is to be noted that the possibility of using a reference state that is uniform is made possible 


by our neglect of gravity in the work from the end of Section III onwards 
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The complete specification of the reference state and the corresponding difference 
is not fixed just by the uniformity requirement but is is a matter of discretion, depending 
on the physical context under consideration. The most natural choice for many purposes 
is to fix the reference state by the requirement that that it give unchanged values for the 
total integrals and averages of the independent current variables, but for some purposes 
other choices might be more convenient. Before adopting any particular convention as to 
the choice of reference state let us first consider the “gross” difference denoted by 6^Q which 
is already well defined a priori. 

In order to see the nature of the dominant contribution to the “gross” differences in the 
large r limit, we exploit the fact that the equation (p.l7|) for the radial dependence of the 
chemical potential will be expressible in the form 




M 


w 


(A12) 


Assuming that the quantity Q{r) under consideration has a dependence on the radial coordi¬ 
nate r only through the intermediary of either /i(r) or r(r), one can use a Taylor expansion 
for Q considered as function of the two variables /i^ and T"^, 


Q(r) = Q(r) + 


dQ 


(p2(r) -/i2(r)) + 


dQ 


(r-"(r) - r-=(f)) + 


(A13) 


where we have written explicitly only the first terms of the Taylor expansion. Using (|3.12|) 


and (|3.15|) , one gets 


M 


II (r) - II (r) = (w - w), 


(AM) 


and 


r-2(r) -r-2(f) 


2Mf22 


/I 

1 \ 

Mc^ 

/ w\ 

\w 

wj 


v-Q 


(A15) 


where the second equality follows from (|5.7|) . 

Taking into account all the orders of the Taylor expansion, one can see easily from the 
previous expressions, that the contribution from any order /3 (/3 > 1) will be a sum of terms 
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of the form w^{w/w)'^ with a < (3. When one substitutes the Taylor expansion (|A13|) of 
Q{r) into the integral in (K3), such a term will yield the contribution of the form 


w 


w'' 


wj a 


a—1 




1 a- 1’ 

except the case a = 1, in which case the contribution will have the form ln{w/xv}. It 
can thus be seen that the expansion for (|A8|) will take the form 


OM = Qir) - !«{-} + o{»} ■ 


(A16) 


whose dominant term is not afected be the rigid rotation whose contribution is proportional 
to and thus at most 0{w}. 

Finally, using once more the Taylor expansion (A13), one arrives at the simpler form 


QW = Q. 


M vldQ 


■uy 


4c^ L/i djjL 

We hnally obtain a limiting relation of the form 

Mw 


w In + 0{w} . 
w 


(A17) 


s^Q~Q 




ln{^}, 

w 


(A18) 


where the “gross” asymptotic deviation coefficient is defined (so as to have the same dimen¬ 
sionality as Q itself) as the logarithmic derivative 

dQ 


Q = 


/i 


dfi 


(A19) 


We now come to the final step, which is to take account of the difference between the 
asymptotic limit value with respect to which the preceeding “gross” deviations are de¬ 
fined, and the corresponding reference state value with respect to which the desired 
“net” deviations are defined. In order to proceed, let us use the abbreviation q to indi¬ 
cates some chosen set (which might consist of the relevant independent longitudinal current 
components) of the independent field variables needed to specify a homogeneous state of 
the system, so that (with implicit sumation over the individual variables in the set q) the 
difference ( |A11|) between the actual asymptotic value and its value in the as yet unspecified 
homogeneous reference state will have an asymptotic form given by an expression of the form 
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AD A S 

'^ooQe ~ ^ ~ ^ (<5oo? - • 


(A20) 


More specificly, let us suppose that the set of variables q is chosen so as to £x a corresponding 
the choice of the reference state in an obviously natural way by simply matching the corre¬ 
sponding average (and therefore also the total) values, so that the homogeneous reference 
state values agree with average values q for the vortex state under consideration, which 
is expressible in terms of the notation of (|A10|) just by 


S^q = 0 . 


(A21) 


Subject to this choice, it evidently follows from (|A20|) that the “net” variation ( |A10|) in 
which we are ultimately interested will be given for a generic quantity Q by 


6^Q^6 Q-^6 q. 


(A22) 


It now follows from the application of (|A18|) to the q variables that there will be a modihed 
net asymptotic deviation coefficient given by 

dQ 


Q = Q - q 


dq 


(A23) 


in terms of which the “net” deviation ( |A24| ) will be given by the asymptotic formula 

^ - Mw 




(A24) 
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